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Abstract. In the paper, the Bismut derivative formula is established for multidimensional 
SDEs driven by additive fractional noise (i < H < 1), and moreover the Harnack inequality is 
given. Through a Lamperti transform, we will show that the Harnack inequality also holds for 
one-dimensional SDEs with multiplicative fractional noise. As applications, the strong Feller 
property is derived and an invariant probability measure is founded for a discrete semigroup. 
Finally, we will establish the Driver integration by parts formulas for SDEs driven by additive 
fractional noise (0 < H < 1) and present the shifted Harnack inequalities. As a consequence, 
the law of solution has a density with respect to the Lebesgue measure. 
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1 Introduction 

The Bismut derivative formula [6] and the Driver integration by parts formula [11] are two 
important tools in stochastic analysis. Let V be the gradient operator and Pj stand for the 
diffusion semigroup. Using the two formulas, one can estimate the commutator VPt — PV, 
which plays a key role in the study of flow properties [14]. Based on martingale method, 
coupling argument or Malliavin calculus, the Bismut derivative formula has been widely studied 
and applied in various fields, for instance, see [12} 133] [36] and references therein. Whereas, by 
using a new coupling argument, [M] established general results on integration by parts formula 
and applied them to various models including degenerate diffusion process, delayed SDEs and 
semi-linear SPDEs. 



On the other hand, dimensional-free Harnack inequality, initiated in [29], has various applica- 
tions, see, for instance, [7] [25] [26] [30] [31] for strong Feller property and contractivity properties; 
[H [2] for short times behaviors of infinite dimensional diffusions; [15] for heat kernel estimates 
and entropy-cost inequalities. In general, one can establish this type of Harnack inequality by 
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using the method of derivative formula (e.g. [U [EU [25l [29]) or the approach of coupling and 
Girsanov transformations (e.g. [H \13 \ |24" 1 [32]) . 

In this article, we are concerned with the stochastic differential equations driven by fractional 
Brownian motion. By using the theory of rough path analysis introduced in [18] . Coutin and 
Qian [8] proved an existence and uniqueness result with Hurst parameter H G (3,5)- Based on a 
fractional integration by parts formula [35], Nualart and Ra§canu [21] established the existence 
and uniqueness result with H > ^. For the regularity results about the law of the solution, one 
can see [16[ \19\ [23] and references therein. Recently, Fan [13] established the Bismut derivative 
formula and Harnack inequality for SDEs driven by fractional Brownian motion with H < A, 
and Saussereau [28] proved transportation inequalities for the law of the solution with H > A. 

The main purpose of our work is to study the Bismut derivative formula and the Driver inte- 
gration by parts formula for SDEs driven by fractional Brownian motion. First, using coupling 
argument and the Girsanov transform for fractional Brownian motion, we will show that, in the 
case of H > ^, the Bismut derivative formula holds for multidimensional equation when the 
diffusion coefficient is constant. Based on the derivative formula, the Harnack inequality will 
be given for multidimensional equation with additive fractional noise, and for one-dimensional 
equation with multiplicative fractional noise. As applications of the Harnack inequality, the 
strong Feller property is derived and an invariant probability measure is founded for a discrete 
semigroup. Finally, we will establish the Driver integration by parts formulas for SDEs driven 
by fractional Brownian motion with H G (0, 1). Consequently, the shifted Harnack inequalities 
and the existence of the density with respect to the Lebesgue measure for solution are presented. 

The paper is organized as follows. In section 2, we give some preliminaries on fractional 
Brownian motion. In section 3, we will prove the Bismut derivative formula for multidimensional 
SDEs with additive fractional noise (H > and moreover obtain the Harnack inequality. 
Through a Lamperti transform, we will show that the Harnack inequality also holds for one- 
dimensional SDEs with multiplicative fractional noise. The remains of the part is devoted to 
establishing the strong Feller property and invariant probability measure. Finally, section 4 
discusses the Driver integration by parts formulas and, as a consequence, shows the shifted 
Harnack inequalities and the existence of the density with respect to the Lebesgue measure for 
solution. 

2 Preliminaries 

In this part, we recall some basic facts about fractional Brownian motion. For more details, one 
can refer to |22j. 

Let B H = {Bf 1 , t G [0, T]} be a ci-dimensional fractional Brownian motion with Hurst param- 
eter H G (0, 1) defined on the probability space (O, J£",P), i.e., B is a centered Gauss process 
with the covariance function KB^' l B^'^ = Rn(t, s)5ij, where 

R H {t,s) = \{t™ + s™-\t-s\™). 

In particular, if H = \,B H is a <i-dimensional Brownian motion. By the above covariance 
function , one can show that ~R\B^' 1 — B^ ,l \ p = C(p)\t — s\ pH , Vj> > 1. As a consequence, B H ' 1 
have (H — e)-order Holder continuous paths for all e > 0, i = 1, • • •, d. 

For each t G [0, T], we denote by Tt the <r-algebra generated by the random variables {Bf : 
s G [0,i]} and the P-null sets. 
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Let <§ be the set of step functions on [0, T] with values in M. d and % be the Hilbert space 
defined as the closure of <§ with respect to the scalar product 

d 

((-f[o,ti]>- ' ''-fa,**])' ( J [o,si]) ■ • -J[o,s d ]))n = y~]RH(ti,Sj). 

i=l 

d R . 

The mapping (L 0ttl j, ■ ■ ■, I[o,t d ]) ^ Yl ' can be extended to an isometry between H and the 

i=l 

Gauss space H\ associated with B H . Denote this isometry by 4> \— > B H ((f>). 

On the other hand, it follows by [TOJ that the covariance kernel Rn(t, s) can be written as 

rtAs 

R H (t,s)= / K H (t,r)K H (s,r)dr, 
Jo 

where Kh is a square integrable kernel given by 

K H (t, s) = T(H + l -y\t - s) H -lF{H - 1 1 - H, H + \, 1 - ~), 
in which F (■,■,-,■) is the Gauss hypergeometric function. 

Define the linear operator K* H :S L 2 ([0, T], R d ) as follows 

(K* H <p)(s) = K H (T, s)<f>(8) + | T (0(r) - <j>( 8 ))^(r, s)dr. 

In particular, (if^ (/[ , tl ], • ■ -,^[0,^])) (0 = (K H {h,-),- ■ ■, K H {t d , ■)). 

By we know that, for all G K^ip) L 2^ 0jT ^ m d^ = (4>,rp)n holds. By the B.L.T. 

theorem, K^j can be extended to an isometry between % and L 2 ([0, T], Hence, due to [3], 
there exists a Brownian motion W such that B H ((p) = W(K^(j)), \/(p G 

According to [TO], the operator K H : L 2 ([0, T], E d ) /^ + ^(L 2 ([0,T],M d )) associated with 
the kernel Kh(-, ■) is defined as follows 

(K H f)(t) := f K H (t,s)f l (s)ds, i = l,--;d, 
Jo 

H+- 

where / + 2 * s the (.ff + ^)-order left fractional Riemann-Liouville integral operator on [0,T]. 
It is an isomorphism and for each / G L 2 ([0, T], R d ), 

(K H f)(s) = I™8*- H ll; H 8 H -h, H<\, 
(K H f){s) = I^-h^s^f, H>\. 

H+- 

Thus, for all h G I 0+ 2 (L 2 ([0, T], M. d )), the inverse operator Kjj is of the following form 

{K H l h){s) = 8 H -^D^ S l- H h', H> 1 -, 
( K H lfl )(s) = s^D^s^D^h, H<\, 
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where Dq + is a-order left-sided Riemann-Liouville derivative, a G (0,1). For more details on 

fractional calculus, one can refer to [27] , 

In particular, if h is absolutely continuous, we get 

(K^h)(s) = s H -hl; H s l i- H h', H<\. 

In the paper, we are interested in the following stochastic differential equations driven by 
fractional Brownian motion on R d : 

dX t = b(X t )dt + <T(X t )dB? : X = x, (2.1) 

where b : R d -> R d and a : R d -> R d x R d . 

Define P 4 /(x) := E/(Xf), t G [0,T], / G ^ ft (R d ), where is the solution to the equation 
()2.ip with the initial value x and ^b(R rf ) denotes the set of all bounded measurable functions 
on R d . For < A < 1, we denote by C A (0,T;R rf ) the space of A-H61der continuous functions 
/ : [0, T) — > R d equipped with the norm 

iii/iiia:= ii/iu + ii/iu, 

where H/IU = sup \f(t)\ and ||/|| A = sup ^feP- 

We will establish the Bismut derivative formula and the integration by parts formula for Pt, 
and moreover present some applications. 



3 Bismut derivative formula 

In this part, we fix ^ < H < 1 and deal with the equation (|2.ip with additive noise. Assume 
the diffusion coefficient a is equal to the unit matrix 1^. 

Next we give some assumptions on the drift coefficient b: 
(HI): b is differentiable, and |V6| < K u |V6(x) - Vb(y)\ < K 2 \x - y\, Vx, y G R d , 
where K\ and K 2 are positive constants. 

The aim of the part is to establish a Bismut derivative formula for Pt which will imply the 
Harnack inequality. For / G ^b(R d ), x,y G R d ,T > 0, we will consider 

„ Pfn r P T f(x + ey)-P T f(x) 

VyP T f(x) := lrni . 

e— >o e 



Theorem 3.1 (Bismut derivative formula) Assume (HI). Then, for each T > 0, / G £$b(R d ), x, y G 
R d , VyPxfix) exists and satisfies 

V y P T f(x) = E (f(X T )N T ) , 



where 



N T = — 

r( 



+(H--) (s H -~> S - ^^[( T -r)V y b(X(r))+y]dr,dW s 

2 Jo \ Jo (s-r)2 +H 



1 \ f T I f s (T — r)V y b(X(r)) - (T - s)V y b(X( S )) 



2 Jo \Jo (s-r)2 +H 
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Proof. We restrict ourselves to the case d = 1 for simplicity. 
The proof will be divided into two steps. 

Step 1. For every e > and y G R, let us introduce the following coupled stochastic differential 
equation 

dX\ = b(X t )dt + dBf - ^ydt, X e Q = x + ey, (3.1) 
dX t = b(X t )dt + dBf , X = x. (3.2) 

Obviously the equation (|3.ip has a unique solution. In view of (|3.ip and (3.2), we conclude 
that XI - X t = ^ey, Mt G [0,T], in particular, Xf, = X T . Let r) t = b(X t ) - b(Xf) - ^y and 
Bf = Bf + Jj 7] s ds Vt G [0, T], then we can reformulate (|3.ip as 

dX\ = &(X t e )cft + dBf. 

Observe that, by the Gronwall lemma, 

Halloo < [(1 + K x T)\x\ + r|6(x)| + ||B H || 00 ]e A ' lT , 

and 

\X t -X s \ < [KiWXWoo + Kx\x\ + \b(x)\]\t - s\ + \Bf -Bf\ 

< d 1 \t-s\ + d 2 \\B H \\ 00 \t-s\ + \Bf -Bf\, (3.3) 

where di = K l e KlT \x\ + (1 + i^iTe^ lT )(Ki|x| + |6(x)|) and d 2 = K ie KlT . 

Next we will show that B H is a fractional Brownian motion under some probability space. 

First note that J Q r] r dr G Iq^ 1 (L 2 ([0, T])). 
In fact, by the definition of rj, we have 

r/ t - Vs = b(X t ) - b(X s ) - {b{XD - b{XD). 

Hence, it follows from the condition of < K\, Vx G R that 

\Vt ~ Vs\ < Kx(\X t - X s \ + \Xt - Xl\) < K x (l\X t - X s \ + ^Ae\y\ 



By \21\ Theorem 2.1], it follows that X have a-order Holder continuous paths for all a G (0,H). 
Therefore, the process rj also have a-order Holder continuous paths for all a G (0, H). It implies 

that n G l"~^{L 2 ([0,T})) and moreover f r] r dr G 1^ (L 2 ([0, T])). 

According to the integral representation of fractional Brownian motion and the derived above 
fact, we deduce that 



Bf = [ K H (t,s)dW s , 
Jo 



where W t = W t + f^K" 1 f Q Vr dr)(s)ds. 
Let 



R e := exp 



£ (k h 1 J Vr dr^J (s)dW s (k- 1 jf Vr dr^ (s)ds 
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Now we are in position to prove that (i?/ f )o<t<T is an ^-fractional Brownian motion with 
Hurst parameter H under the new probability R e P. Due to the Girsanov theorem for the 
fractional Brownian motion (see e.g. |10l Theorem 4.9] or [201 Theorem 2]), it suffices to show 
that ER t = 1. 

2' yli c -'OH 



Recalling the facts: in case of H > |, Mh G I 0+ 2 (L 2 ([0, T])), 



{K^h){s)=s H - 1 >D^{s*- H h>), 
and V/ G iJ+fLP), < a < 1, p > 1, 



»+'' r(i - a) V i° Jo (i - </) 1+ " 



we can get 



(k- 1 J^ r dr\s) = s H -^D^{r^- H Vr ){s) 



r(|- 




i 




r(l- 




+(#- 


-i) 


i 




r(|- 





1 (.i-^ +( H-v-t r sh " ,h \z H H ^ A 

\ 2 Jo (s-r)2+ H J 

S2 n rj s + {H )s 2 / rj r dr 

2 Jo (s-r)2 +H 

dr) 

o (s-r)2 +H J 
(J 1 + / 2 + J 3 ). (3.4) 



Notice that \rjt\ < (K\ + A)|y|e, and 



dr= 11 — _X.de ■ s l ~ 2H =-.c s l - 2H 



o (s-r)2 +H Jo (1 



+H 



where we use the change of variables 9 = - for the first integral and Cq is a positive constant. 
What remains to be dealt with is the term I3. Let F(x,y) := b{x) — b(y), we have 

Vs ~ Vr _ dr _ f s F(X* r ,X* s )-F(X r ,X s ) dr 



(s — r)i +H Jo (s — r)^ +H 



^FjJ^JQjXl - X r ) + d 2 F(X r ,X s )(XI - X,) ^ 

(s — r)z +H 

e f s b'(X r )(T-r)-b'(X s )(T-s), 
= —y / j dr 

T Jo (s-r)2 +H 

e f s (b'(X r ) - b'(X s )){T -r)+ b'(X s )(s - r) 

1 Jo (s-r)2 +H 
where we use the mean value theorem for the function F in the second relation and X r 

x r + o^xt - X r ), x s = x s + o x {xt - X s ), < Q x < 1. 

By (HI) and (|3.3p . we obtain 



[' dr 
'0 (s-r)2 +H 



< K 2 \y\e [ S } Xs ^[ dr + KMe !\s-r)^ H dr 
Jo (s-r)2 +H J- Jo 
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< K 2 \y\e 



fe(s-r) + \X s -X r \ 



(s — r) 



(s — r)' 



-u 



dr 



o 



< 



rdiif 2 + ^2|y|e + ^i 



+ 1 ^| y | a! H||J3*|| Jf _ 
2 " 

where 8 is taken such that < 5 < J. 

Now, we can estimate Kjj 1 (J rj r dr)(s) as follows 



if) T 



|y|es 



^- H + ^%\yM\B H \\oosfi- E 

2~ H 



Krr 1 / r/ r cir 



< des^ H + C 2 e S i- H + C^es^WB 11 ^ + C^a* - *!! 5 * Ilff- 



where 



Ci 



m-H) 



l + C [H-- 



c 2 



1 H -^Td 1 K 2 + K 2 \y\e + K 1 



\y\, 



c 3 = 
c 4 



1 



r(|-^)3_ H 
i if- 



H — - 

2 d 2 K 2 \y\, 



\-5 



K 2 \y\- 



Moreover, the above estimate follows that 



o 



< Ae z 



r\ r dr I (s) 



o 



-2H 



ds 



+ 



2-2H 4-2H 
As a consequence, we obtain 



+ 



°4 J || dJ? ||2 

2-25 ^ ^ 



E exp 



< C 5 Eexp[e 2 (C 6 \\B H \\l c + C 7 \\B H \\ 2 H _ 

< C 5 Eexp[(C 6 VC 7 )e 2 \\\B H \\\ 2 H _ s ] 



where C5 = exp 



2e 



2 f C 2 T 2 ~ 2H C£T 



2-2H 



4-2H 



c fi 



and C 7 



CfT 2 



(3.5) 



(3.6) 



2-ff w ' ~~ 1-5 • 

Using the Fernique theorem, we know that, the right-hand side of (I3.6P is finite when e is small 
enough. Consequently, the Novikov condition yields that KR e = 1. 

Step 2. From step 1, we have already known that, for small enough e, {B^ )o<t<T is an 
fractional Brownian motion with Hurst parameter H under the probability R e P. So, under R e P, 
the law of the process (X^' x ey )te[o,T] is the same as 

te[o,T] under P, where the superscript 
x + ey stands for starting part of the corresponding equation. By the fact: X^ x+ty = X^, we 
get 



P T f(x + ey) = Ef(X* +ey ) = E(RJ(X^ x+ey )) = E(i? e /(Xf )). 



(3.7) 
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Hence, (|3.7p implies that 



V,P T / W : =1 . m F ^ + ^- Pr/W = l im 



E/(X X ^ e 



iL - 1 



(3.8) 



Define M t := - jj (ff^ 1 / ?? r dr) (s)dW s , then i? e = exp[M T - |(M) T ]. Next we will show that 



limE I /(Xf)^^ ) = limE ( /(Xf) 



(3.9) 



Since / is bounded, it suffices to prove that 

limE gi^ = limE ^-^. 

Without lost of generality, in the remaining proof, we suppose that e < 1. 

Indeed, noting that, for all x 6 R, |e x — 1 — x| < x 2 e' :r ' and x 2 < e^' hold, we deduce that 



Re - 1 - (M T - UM) t ) 



1 



1 



< - M T --(M) T e 



\M T \ + ±{M) q 



i / 1 

e5 ( ^M T 



< e h^ lMTl+ ^ {M)T e^ + ^ 

i / ^ Mt --^(M) t +e ^M T -J^{M) T ^ c -^(M) T +-^(M) T 



Using the Holder inequality, the C r inequality and (|3.5p . the above inequality yields that 
Re - I - (M T - l(M) T ) 



E 



: el f 2 Ee^ MT -^ (M>T 



+ 2Ee 



(4e)5 (Ee^ 



J^( M )T + Jl2 (M)t\ ^ 



which shows that (|3.9p is true. 

Noting that ([33]) and Edl^H^ + ||-B J3 ]|l r _ < y) < oo, we may write as follows 



limE ( /(Xf) — = limE [ /(Xf) — 

e->o V e / £->o V e 



(3.10) 



In view of (|3.4p . Mt can be decomposable as 



T (k- 1 J Vr dr ) (s)dW ; 

i r 



z JO 



(s - r) 



rj r drdW s 
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T rs 



+(-l)(H--) 



o Jo (s - r) 



2—drdW a 
i i it * 



r(|-fl) 



[ h + J2 + J3 



(3.11) 



For the term J\ , making use of the Jensen inequality, we obtain 



E 



< 







s \-»(l + b'(X s ){T-s))V-dW s 







, T = i-2H ( KXl) - b(X s ) - bVQ^ey 



ds 



E l T s l - 2H {b'{X s +e 2 {Xt-X s ))-b'{X s )) 2 (^— -y) ds 



(3.12) 



where we use the mean value theorem in the last expression and #2 £ (0, 1). 

The dominated convergence theorem implies that the right-hand side of (|3.12p goes to as e 

tends to 0. 

Similar to the treatment of Ji, for J2, we conclude that 



E 



T 



1_I7 1 — Tt 

1 / " s2 n — ra 11 



s 2 



(s — r) 



--±-(l + b '(X r )(T-r))^ 



drdW, 



0, (3.13) 



as e tends to 0. 

For the term J3, we get 



E 



T r s 



r] s -rj r 1 



T rs u 



J (s — r) 



&'(X r )(T-r)-6'pQ(T- S )y 



Jo 



(s — r) 



drdW, 



< < E 



E 



b(X s ) - b(X r ) - [b(Xp - b(X< r )] - [b'(X r )(T - r) - b'(X s )(T - S )]f ^ 

e(s-r) l 2+ H 

[b'{X r ) - b'(X r ) - b'{X s ) + b'(X s )](T - r) + [b'(X s ) - b'(X s )](s - r) y 



ds 



(s - r) 



T 



dr 



ds 



(3.14) 



where we apply the mean value theorem in the last relation and X r = X r + 9^{X^, — X r ), X s 

x s + e 3 (xt-x s ), 3 e (0,1). 

We first claim that 



lim 



\b'{X r ) - b'(X r ) - b'(X s ) + b'(X s ))(T -r) + [b'(X s ) - b'(X s ))(s - r] 



(s - r )2 +H 

In fact, noting that X s — X r = X s — X r + 6*3 -^ey, by (|3.3p we get 
\b'(X r )-b'(X r )-b'(X s ) + b'(X s )\ 



dr = 0. (3.15) 



(s — r) 
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< K 2 



< K 2 



< K 2 



\X S — X T \ + \X S — X r 



(s-r) 



2\X s -X r \ + s -^\y\ 



(s-r) 



2d 1 (s - r) + 2d 2 \\B H \\ O0 (s - r) + 2\B? - B?\ + Sf\y\ 



(s — r)- 



-ii 



s-r)f- H + 2i^ 2 C(o;)(s 



-J3 



where we use the fact: 2?^ have /3-order Holder continuous trajectories for all < j3 < H and 
choose \ < (3 < H here. 

Consequently, (|3.15p follows from the dominated convergence theorem. 

In order to prove (|3.14p converges to 0, as e tends to 0, combined with (j3. 15|) . it suffices to give 
a majorizing function due to the dominated convergence theorem. 
Notice that 

[b'(X r ) - b'(X r ) - b'(X s ) + b'(X s )](T - r) + [b'(X s ) - b'(X s )](s - r) y 



(s-r) 



-dr 



< W 



2di + M + 2d a ||J3- H 'I|oo ) (s - r)*- H - 



r>H r>H\ 
i>„ — D„ 



(s - r) 



T 

2 



dr 



2#i | 



(s — r) 2 ^cir 



^ ^ „3-2H , n e 3-2iJ || R i?||2 j_ r 

< G 8 s + o 9 s LB Hoo + ^ios 



1-2(5 1| o-ff 1 1 2 

\\B \\h-S> 



where Cg, Cg and Cio are positive constants. 

Due to flSHZD , ([3II5D and (ETT41 . the proof is complete. 

Remark 3.2 The Bismut derivative formula presented in Theorem 13.11 can be easily extended 
to the following equation 

dX t = b(X t )dt + cr(t)dBf, X = x, t G [0, T], 

where the stochastic integral exists pathwise under proper assumptions. 



As applications of the Bismut derivative formula derived above, we may get the explicit gra- 
dient estimate and the dimensional free Harnack inequality for Pj-- 

By the Fernique theorem, there exists a positive constant Ao such that Bq := Eexp[Ao|| \B H \ \\%_g] < 

oo. We set ||VP T /(x)|| := sup \V y P T f(x)\, Vx G R d . 

y m d ,\y\<i 

Corollary 3.3 Under assumptions of Theorem \3.1[ Let f G ^^{W 1 ) be positive. Put 



A 



l •= 



(T(§ - H)fT™ 



(1 + KiT) 2 + [C (22-i)(l + 2f 1 r)] 2 + 



2H 



2H 



V3(H- \)(d x K 2 T + K^T 
V4 - 22/(1 - H) 



A 2 :-- 



V(H-W 
(T(l-H)y 



d 2 K 2 



^-2H 



2H 



V 



T 



2-2(5 



25 
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(1 ) For each x£K d and a> J , we have 



||ViW(aO|| < a[P T (f log f)(x) - (P T f)(x)(logP T f)(x)] + - ( A 1 + ^logB 



(2) Let p>l. Then, for all x, y G R d such that \x - y\ < ^y^, we get 

A 



(P T f(x)) p < P T f p (y)exp 



^(A 1 + ^logB )\ X - y ? 



(3) For each x G R d , lim Prf(y) = Prf(x), i.e. Pr is strong Feller. 

Proof Let / G &b(M. d ) be positive. Combining Theorem 13.11 with the Young inequality (see 
e.g. [5j Lemma 2.4]) yield that, for each a > 0, 

|V„Pr/(x)| < a[P T (f log f)(x) - (P T f ) (x) (log P T f)(x)} + [alogEe^]P T f(x). (3.16) 

Now we turn to calculate Eea*. 
By the expression of Nt, we obtain 



(N) T < 



[r(| - h)T\* 



A-2H\ 



-'\(T-s)V y b(X(s)) + y\ z ds 



1. 



T 



2H~l 



e 1 TT 1 IT 



o [s — r p 



[(T-r)V y b(X{r))+y]dr 



ds 



■ \T(\ - H)T] 2 
A simple calculus shows that 



T 



(T - r)V y b(X(r)) — (T — s)V y b(X(s)) 



(s — r) 



dr 



ds 



{h + h+h}. 



(3.17) 



h < 
h < 

h < 



2-2H 



-T^\y\\ 



[C Q {H-I){l + K X T)1_ T2 _ 2H ^^ 



2-2H 



A-2H 
3 

+ 2-25 



(H - \){dxK 2 T + K x ) 



\-8 



^4-2H|„,|2 



\y\ z + 



4-2H 



(H - \)d 2 K 2 
l-H 



T^l^Wllyf 



T 4 - 2S \\B H f H _ s \y\\ 



Therefore, we get 



{N ) T <( Al +A 2 \\\B H \\\ 2 H _ 
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I n A 

Then, for any a > w =: ckq, we have 











Eexp 




< ^Eexp 






a 




or 



< exp 
= exp 
= exp 



Ai, |2 

—\y\ 

cr 

— \y\ 



Eexp 



2^2 ,.,2 



7/HII R^lll 2 

y\ \\\^ \\\h-5 



4 

2a 2 



a 



Eexp [A |||5 H ||&_ 5 ])' : 

21 



A 1 + ^logB ^ 



Ar 



(3.18) 



(1) By (I3TT6D and (|Q5j> . we get, for a > y^|f|y 



|V y P T /(x)| < a[P T (/log/)(x) - (P T /)(x)(logP T /)(x)] + - Ai + log Bo |yrW(x) 



a 



A, 



A 



Let \y\ < 1, we derive the desired result. 

(2) The part of the proof follows the arguments of [5j Theorem 1.2]. Let f3(s) = l+s{p—l), r y(s) 
x + s(y — x), s S [0, 1]. We obtain 



m p Tfm(l{s)) 



(7W) 

x 



p — 1 

^7 



Pt(/^ (s) log/^ fl >) - (Pt/ (s) ) logPr/^l (7(a)) " V^Pt/^tOO) 



Then, for x, y € M d satisfying \x — y\ < 2 — J ^P-, it follows from (|3.16|) and (|3.18|) that 



p \ 2A 2 - 



> 



p 



fi{s)a 

P 
p-1 



Ai + ^logPo) \y-x\ 2 

A 2 

Ai + — logP ] \y- x\ 



where we take a = ^t4 in (|3. 16[) . 

Integrating on the interval [0, 1] with respect to ds, we get the desired result. 

(3) According to the proof of [9] Proposition 4.1], the above result (2) implies the strong Feller 

property of Pp. 

Remark 3.4 The Harnack inequality in Corollary \3.3\ is local in the sense that \x—y\ is bounded 
by a constant. How to establish a global Harnack inequality is an interesting problem. 

The Harnack inequality for one-dimensional equations with multiplicative noises can be de- 
rived from Corollary 13.31 Now we assume that d = 1 and consider the stochastic differential 
equation 



dX t = b(X t )dt + u{X t )dB?, X = x, t€ [0, T]. 
We give the following assumptions on the coefficients: (H2) 



(3.19) 
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(i) b and a are both twice differentiable, and the functions themselves and their derivatives 
are bounded; 

(ii) There exist two constants d% and d± such that d^ > d^ > and d$ < a(x) < d^, \/x G R. 

Under the above assumptions, |21j shows that there exists a unique adapted solution to equation 
(|3.19p whose paths are Holder continuous of order H — e for every e > 0. Note that the stochastic 
integral appears in equation f|3. 19j) can be considered as pathwise integral. See [35] for more 
details. Set Pt/(x) := E/(X^), where (Xf )t^[o,T] is the solution to equation (13.191) with initial 
value x. 



Theorem 3.5 Assume (H2) holds. Let p > 1 and f G ^(R) be positive function. Then, there 
exists two positive constants Ai and A 2 such that the Harnack inequality 



{P T f{x)Y < Pr/ P (y)exp 



holds for all x,y G R satisfying \x — y\ < T -^rd^ x / tta~ ■ 

Proof Put F(y) := J Q V -^^dz. Then, according to the change-of-variables formula [351 Theo- 
rem 4.3.1], we know that X is the unique solution to equation (13. 19ft if and only if the process 
Y = F(X) is the unique solution of 



dY t 



b(F- 1 (Y t )) 



dt + dB t H , Y = F(x), t G [0,T]. 



(3.20) 



Define Prg(z) := Eg(Y^), where (^ 2 )fg[o,T] is the solution of equation (|3.20p with initial value 
z. 

Note that -^F~ l {x) = a(F~ 1 (x)). Hence, due to the assumption, we deduce that the coefficient 
fe l^-il X NN is twice differentiable and moreover, 



b^F-^W b'a-ba' 



<F^x))J a ^ 

bjF-ijx)) }" _ b"a 2 - baa" - bW + ba' 2 
°{F-^x))) ~ a 



(F-\x)) 



are both bounded. Then, by (2) of Corollary 13.31 we conclude that there exist two positive 
constants A\ and A 2 such that, for each g G <5^,(R), 



{P T g{zi)) p < P T g p (z 2 )exp 



p 



p — 1 



Ao 

A! + — log B ) \z x 
^0 



z 2\ 



(3.21) 



holds for all z\,z 2 G R satisfying \z\ — z 2 \ < ^y- J 7^ ■ 

For x,y G R such that \x - y\ < ^d 3 ,/^ and / G ^&(R), we take z\ = F(x),z 2 = F(y) 

and g = f o Then applying (|3.2ip to the particular choices z\, z 2 and g, we get the desired 
result. Indeed, first we have 



\zi - z 2 \ = \F(x) - F(y)\ < - y\ < - — -J -2- 

«3 P V 2A 2 
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Furthermore note that Y t F{x) = F(Xf), t G [0,T]. Then, we obtain 

Prg( Zl ) = M 9 (Yt) = E / ° F-\F(X£)) = E/(Xf ) = P T f(x) 

and 

P T ^(z 2 ) = EgP(F(X^)) = E[f o F-\F(X^))f = Ef{X y T ) = P T f p (y), 
which show that the proof is complete. 

Note that the solution X of equation (|2.ip is not a Markov process. Consequently, (Pt)t>o 
does not consist of a semigroup. So, we introduce the following semigroup in discrete time, i.e. 
for any Borel set A in M. d , 

P T {x,A):=P T I A (x), P$(x,A) := f P^Hx, dy)P T (y, A), n > 2. 
In general, (P$f)(x) = f Rd f(y)P$(x,dy), x G M. d , f G ^ 6 (]R d ). 

Next let g = 1^ in the equation ()2. 1|) and we will concern with the existence of invariant 
probability measure for semigroup (-Pp) n >i for fixed T > 0. Besides the assumption (HI), we 
will introduce another important condition. 

(H3): (x,b(x)} < K 3 \x\ 2 , where K 3 £ R satisfies that C 12 e 2K ' iT < 1, C l2 is a positive constant 
given in the proof below. 

Theorem 3.6 Assume (HI) and (H3). Then, the semigroup (P^) n >i has an invariant proba- 
bility measure. 

Proof We will make use of Krylov-Bogoliubov's method. 
Let xq G M. d and define 



it 



Y $x Pt 



k=l 

\i n := , n > 1, 

n 



E P£f(xo) 



i.e. for each / G & b (R d ), /J, n (f) 
Next we will prove the tightness of {/i n }n>i- 

By the change-of- variables formula |35} Theorem 4.3.1] and (H3), we get 

fT rT 



X X T \ 2 = \x\ 2 + 2 f (Xf,b(Xn)dt + 2 [ (X?,dBf) 
Jo Jo 

<\x\ 2 + 2K 3 \X?\ 2 dt + 2 (X?,dB?). (3.22) 
Jo Jo 



'0 Jo 

According to |28t Lemma 5], we get the estimation of pathwise integral in (|3.22p as follows 



o 



T T 

< -^i\\B H h (ll^lloo^ + \\X\\ P T^) , (3.23) 



where L(> 0) and j3 G (|, H) are two constants. 
By (|3.3p . we easily get 



l*IU < \x\ 2 + ( (1 + ^p 2e2XlT + |&(*)|Te^) + e^ll^lU, 
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\\X\\p < drT 1 ^ + d 2 T 1 ^\\B H \\ 00 + \\B H \\ P . (3.24) 
Substituting (I3.23P and (|3.24p into (I3.22|) and taking expectation give 



E|Xf| 2 < Cu + C l2 \x\ 2 + 2K 3 / E|Xf' 2 

Jo 



dt.. 



where 



C 



11 



2LT@ 



d x T + (1 + K ^) 2 e 2KlT + |6( X )| T ^ nB Hy 



+(e KlT + ^Edl^lUH^I^) + T?E\\B H ||| 
Ci 2 = 1 + ^EH^II^. 

The Gronwall lemma yields 

E|Xf| 2 < (Cn +Ci 2 |x| 2 )e 2 ^ 3T =: Ci 3 + Ci 4 |x| 2 . (3.25) 

Now we will show that {J* K d |a;| 2 P^(xo, dx)} is bounded by using an induction argument. 
When 77, = 1, it follows from (|3.25p that f Rd \x\ 2 Pr(x$, dx) < C\ 3 + Ci4|x| 2 . Suppose that 

f ixfp^ixo, dx) < Ci 3 (1 + C14 + c 2 4 + • • • + c^ 2 ) + q^M 2 

JR d 

holds, then by (|3.25p we have 

f \x\ 2 P$(x ,dx) = [ [ \x\ 2 P^ 1 (x ,dy)P T (y,dx)= [ E^P^^dy) 

JR d JR d JR d JR d 

< C 13 + C 1A f d \y\ 2 P^-\x Q ,dy) 

< C13 + C14 [C13 (1 + C14 + C 2 4 + • • • + C™- 2 ) + Cj^Vol 2 ] 

= C13 (1 + C 14 + C 2 4 + • • • + C™- 1 ) + C^lxor- (3-26) 



Therefore, we have, for any n > 1, 



\x\ 2 P?(x ,dx) <-^ + \ Xo \i, (3.27) 

1 - U14 



and moreover, 

n 

E / H d |x| 2 P|(x ,dx) 



|z| 2 /i„((fx) = — < — ^ h l^ol 2 - 

n 1 — C14 

Using the Chebyshev inequality, we obtain 

1 ( C 3 \ 

sup/i n (| • | 2 > r) < - I X — + |x | 2 ) ->■ 0, r ->■ 00, 

n r \ 1 — U14 / 

which shows the tightness of {/Li n }n>i- 

So, by the Prohorov theorem, there exists a probability /i and a subsequence fj, nk such that 
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fJ'nk ~^ A 4 weakly as — > oo. For simplicity of notation, we denote \i n — > /i weakly as n — > oo. 

Now we will prove that /j is invariant for (-Rp) n >i- 

Denote C&(R ) by the set of all bounded continuous functions on M^. 

For any / G C fe (IR d ), by (3) of Corollary O we know that P^/ G C b (M d ), Vra > 1. Moreover, 
we deduce that, for all / G N, 

Hiftf) = lim fi n (Py) 

n— >oo 

n 

E Pr +l f{^) 

v fc=l 

= lim 

n n+l I 

E J?7(so) E ^F7(*o) E JT/(a:o) 

T 771=1 , T 771=71+1 771=1 

= Inn h lim lim 



71 71— >oo n n— >oo 77, 



The proof is complete. 



4 Integration by parts formula 

In this part, we will establish the Driver integration by parts formulas for SDEs driven by 
fractional Brownian motion. As a consequence, the shifted Harnack inequalities are presented. 
In contrast to known Harnack inequality, in the shifted Harnack inequality, a reference function 
is shifted rather than the initial point. 

Now consider the following SDE with additive fractional noise on W 1 : 

dX t = b(t, X t )dt + dBf, X = x. (4.1) 
Throughout the section, we will give some regularity assumptions on the coefficients &(•, •) : 

[o, t] x m d -> m d . 

(i) If H < \, the mapping t ->■ b(t,0) is bounded on [0,T] and \Vb(t,-)(x)\ < K x , Mt G 
[0,T], x G R d , 

(ii) If H > 2i ^ is independent of the time variable, and b is differentiable on space variable, 
|Vb| < K u \Vb(x) - Vb(y)\ < K 2 \x - y\, 

where K\ is positive constant and K±, K 2 are the same as (HI). 

By [20] and [21], we know that the equation (|4.ip has a unique solution. 

Theorem 4.1 Let y G M d be fixed. For all f satisfying f,V y f G & b (M. d ), we get, 

(1) ifH < §, P T (V y f)(x) = E{f{X?)N^), where the random variable is given by 

N r = fjp-jr f (•""* 1 ^J^ 1 ' 9 " rV »' (r ' ■ )(X(r))ldr ' ; 
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(2) if H > \, PT{Vyf){x) = E(f(X£ )N T ), where the random variable N T is shown by 



Nj, = — -5- 



+ (#--)/ Is 11 '* ^^[ y -rV y b(X(r))]dr,dW s 

* Jo \ Jo (s — r)2 +n 



'0 (s — r) '■■ 
7 yb{X{r)) - 
2 Jo \Jo (s - r)^ +H 

Proof Without lost of generality, we show the result for d = 1 for simplicity. 
For fixed initial point x, let Xt solve (I4.ip . On the other hand, for each e G [0, 1], let X\ solve 
the following equation 

dX\ = b{t, X t )dt + dB? + ^ydt, X e = x, te [0, T]. (4.2) 
It is clear that X\ = Xt + ^y, t G [0, T]. In particular, = Xt + ey. 

Define & := &(t,X 4 ) - 6(i,X|) + and £?f = + / * £ s <is. Then we may rewrite ([42]) as 
follows 

dX t e = 6(t, X t e )(ft + dSf, X% = x, te [0, T]. (4.3) 
Note that J £ r dr G (L 2 ([0, T])). 

Indeed, if -ff < |, it suffices to show that J* Q T ^dr < 00, a.s. If H > |, a similar argument to 
step 1 of Theorem 13.11 can deduce the desired result. 

Therefore, according to the integral representation for fractional Brownian motion, we obtain 



B? = f K H (t,s)dW s , 
Jo 



where W t = W t + J^K^ 1 j Q ^ r dr)(s)ds. 
Now let 



R e := exp 



£ (k^ J Zrdr^j {s)dW s -l£ (k^ J £ r dr\ {s)ds 



Following the argument of |13} Theorem 4.1] (H < ^) and Theorem 13.11 (H > ^), we may 
conclude that the process B H is a fractional Brownian motion under the probability R e P, i.e. 
(X, X e ) is a coupling by change of measure with changed probability R t P, and moreover the 
relations ^ Id 1 

— R e \e=o = —N T , H < -; — R e \e=o = —N T , H > -, 

holds in V-(P). 

Hence, due to [MJ Theorem 2.1], we get the desired results. 

As the Bismut derivative formula implies the Harnack inequality, we also deduce the following 
shift Harnack inequality from the Driver integration by parts formula. 

Corollary 4.2 Under the assumptions of Theorem \4.1[ For any f G ), we have the 

following result: 
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(1) ifH<\, then, for each y G R d , 

(p T /) p <(M/(y + -)F)ex P 



p ( B (l-H,±-H)\ (1 + JfxT) 2 



p-iI r(i-ir) 



T 2H 2{l-Hy 



(2) if H > |, then, for every y G K d satisfying \y\ < ^y-J^, 



(PTf) p <(P T {f(y + -)Y)^P 



p 



A 1 + ^logB )\y> 2 



Proof Note that 



(T( 



r 



< 



|y| 2 / s 2 "" 1 



r 2 



o (s — r)2 



1 + KiT 

r(|-PT)T / 

P(f-P,l-ff) V (1 + i^T) 2 2 
r(i-F) i T 2 "2(l-tf) ' 



--H "1 

^— Y— [y - rV y b(r, -)(X(r))]dr I ds 
o (s-r)i + J 

2 



dr 



ds 



and 



(iV 2 )r< (^i +^111^111^)12/ 



where A\ and A2 are defined before. By [341 Proposition 2.3] and Theorem 14.11 we get the 
desired result. 

An important application is to show the existence of the density with respect to the Lebesgue 
measure for solutions of (14.11). 



Corollary 4.3 Under the assumptions of Theorem \4.1[ If H < ^, then the law of the random 
variable Xj, has a density with respect to the Lebesgue measure on H d . 



Proof Let 



C 
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p B{l-H,\-H)\ (l + K^f 



p-i\ r(i-ij) 



J T 2H 2(1-H)' 

By (1) of Corollary S3 we get, for any / G £B b (HL d ), 

{P T f{x)fe~ c ^y\ 2 < (P T {f(y + ) (x) 

Taking / = I a, A is Lebesgue- null set on and integrating both sides with respect to dy on 
the above relation, we deduce that 



(P T I A (x)) P [ e~ c ^y\ 2 dy < 0, 



which implies the desired result. 



18 



Remark 4.4 In \34\ . the author showed that the study of the Driver integration by parts 
formula and shift Harnack inequality was in general more difficult than that of the Bismut 
derivative formula and Harnack inequality. Moreover, he gave some applications of integration 
by parts formula and shift Harnack inequality, for instance, to estimate the density with respect 
to the Lebesgue measure for distributions and Markov operators. 
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